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IMPLEMENTATION  OF  A  REAL-TIME  DATA  PROCESSING  SYSTEM 
FOR  RADAR  BANDWIDTH  EXTRAPOLATION  AND  STRETCH 


Introduction 

A  program  is  underway  at  NRL  to  demonstrate  the  utility  of  wideband  waveforms  in 
radar  surveillance.  Current  effort  has  focused  on  the  development  of  a  demonstration  system 
that  implements  long  range  target  recognition.  The  system  uses  NRL’s  Senrad  radari  which 
performs  normal  air  surveillance.  On  targets  of  interest,  a  wideband  waveform  replaces  a  few 
pulses  of  the  normal  dwell.  Since  the  range  of  interest  is  limited,  these  pulses  are  processed 
using  a  “stretch”  technique  to  yield  high  resolution  range  profiles  of  the  target  Combining 
range  profiles  with  tracking  data  makes  it  possible  to  classify  uncooperative  targets. 

Usually,  increased  bandwidth  will  improve  the  performance  of  target  classification. 
However,  in  practical  radar  systems,  the  amount  of  av^lable  signal  bandwHh  is  limited.  It  has 
been  suggested  that  bandwidth  extrapolation  processing  can  provide  eome  of  the  benefits  of  a 
greater  bandwidth  without  the  problems  of  actually  implementing  the  waveforms.2  Therefore, 
the  objective  of  this  project  was  to  add  real-time  processing  for  bandwidth  extrapolation  to  a 
stretch  data  processing  system. 

Background 

It  is  necessary  to  explain  the  operation  of  the  existing  system  in  order  to  show  the 
necessity  of  this  change  and  to  explain  the  requirements.  The  existing  system  captures  radar 
data  from  the  Senrad  receiver  and  does  the  processing  shown  in  Rg.  1  using  a  general  purpose 
Hewlett  Packard  (HP)  computer.  The  scan  rate  of  the 


N  Pulses  10CX)  Samples  Range  Profiles 


Rg.  1  -  Block  diagram  of  stretch  data  processing 


antenna  is  approximately  five  secorKls  and  the  processing  was  just  completed  before  data  from 
the  next  scan  was  available.  Since  bandwidth  extrapolation  requires  substantial  computation, 
the  system  is  incapable  of  running  in  real-time  without  some  modifications.  The  modification 
deckM  on  was  to  have  some  of  the  processing  done  on  a  different  processor  which  would 
interface  with  the  existing  system.  The  new  system  would  capture  live  data  from  the  radar  also 
and  perform  the  bandwidth  extrapolation  and  stretch  processing  on  this  data.  The  results  would 
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be  sent  to  the  original  HP  computer  where  the  processing  would  be  finished  and  the  results 
would  be  displayed.  The  resulting  processing  times  of  both  the  new  processor  and  the  HP  would 
allow  the  process  to  run  in  real-time. 

The  changes  necessary  to  meet  this  objective  included  translating  the  previous  code 
written  in  the  HP  Basic  programming  language  by  George  Unde  and  his  group  to  the  C 
programming  language.  The  processing  was  d^  on  a  TMS320C30  (C30)  digital  signal 
processor  (DSP)  1C.  The  introduction  of  the  C30  was  most  easily  achieved  by  interfacing  with 
the  hardware  of  the  Volume  Surveillance  project  This  prqect  is  still  in  progress  and  there  is 
no  published  literature  on  the  project  to  date.  The  processors  used  in  this  system  are  the  same 
as  those  used  in  the  Point  Defense  systam.M«s 

The  existing  hardware  system  used  to  do  the  stretch  processing  consisted  of  a  NOVA 
computer,  an  HP  computer  and  other  hardware  which  wiD  not  be  discussed  here.  The  NOVA 
received  digitized  data  from  the  Senrad  receiver  and  sent  it  to  the  HP  for  processing.  The  HP 
performed  the  stretoh  processing  on  a  set  number  of  pulses,  as  well  as  some  other  processing, 
and  displayed  the  results.  The  information  cfisplayed  included  the  range,  speed,  asp^  angle, 
length,  radar  cross  section  and  identification  of  the  target,  the  signal  magnitudes  of  two  pulses, 
range  profiles  for  each  pulse,  an  average  range  profile,  the  peak  level  and  range  offset  from  the 
profiles  and  the  places  where  propeller  modulation  was  found. 

Referring  to  Rg.  1,  a  set  number  of  pulses  consisting  of  1000  samples  each  is  digitized. 

A  complex  multiply  between  the  first  and  last  pulse  of  the  data  and  an  FFT  are  completed.  The 
range  offset  between  these  pulses  is  then  determined.  The  digitized  data  is  also  amplitude 
weighted  and  phase  compensated.  An  FFT  is  taken  and  the  range  offset  found  previously  is 
applied  to  each  pulse  to  time  align  all  of  the  pulses.  To  complete  the  processing,  the  magnitude 
average,  peaks  and  any  propeller  modulation  in  tire  range  profiles  are  found.  The  logarithm  of 
the  range  profiles  is  taken  and  displayed. 

Hardware 

The  hardware  for  this  project  consisted  basicaiiy  of  three  parts:  the  Volume 
Surveillance  system,  the  interfece  and  the  existing  HP  system.  The  Volume  Surveillance 
hardware  included  a  Sun  computer,  the  SP  (signal  processing)  boards  and  the  A/D  (analog  to 
digital)  converters.  The  existing  HP  hardware  mduded  the  HP  computer.  NOVA  computer  and 
other  components  which  were  needed  in  the  existing  system  but  will  not  be  discussed  here.  Each 
of  the  parts  will  be  discussed  in  turn. 

In  the  Volume  Surveillance  system,  the  SP  boards  (only  one  of  which  is  used)  each 
contain  two  SP  nodes  and  a  transputer.  Each  SP  node  contains  one  C30.  The  transputer  has  four 
bi-directional  communication  links  and  is  responsible  for  routing  the  data  from  and  to  the  C30s 
on  the  SP  board.  A  separate  network  of  transputers  is  configured  to  route  data  between  the  Sun 
computer  and  the  SP  boards.  The  block  dagram  of  the  modified  Volume  Surveillance  hardware 
is  shown  in  Fig.  #2.  The  blocks  marked  SR  (server)  are  the  available  SP  boards.  Only  one  SP 
board  was  used  since  only  one  C30  was  needed  for  processing.  The  blocks  marked  RT  (router) 
and  CS  (console)  are  part  of  the  network  of  routing  transputers. 
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Rg.  2  <  Block  diagram  of  modified  Volume  SurveiHance  hardware 


Some  other  modifications  to  the  Volume  Surveillance  hardware  were  necessary  in  order 
to  capture  the  digitized  data.  The  program  of  the  XiHnx  field  programmable  gate  array  (FPQA) 
1C  on  the  one  SP  board  used  was  modfiied  by  Mike  Livingston.  This  1C  expected  an  external 
triggar  for  each  pulse  when  the  radar  data  was  vedid.  When  one  of  these  triggers  was  rec^ved, 
the  Xilinx  1C  collected  and  stored  1024  digitized  samples  from  the  A/D  converters.  When  the 
data  frx  all  the  pulses  was  ready,  the  software  cfiscussed  below  was  able  to  (firecdy  access  the 
data. 


The  second  modification  was  a  new  A/D  intertace  board  which  was  also  completed  by  Mite 
Livingston.  This  board  buffered  the  digitized  I  and  Q  channel  data  from  the  radar  and  converted 
it  to  a  32  bit  data  format  The  hsvdwwe  used  to  reset  the  SP  boards  ¥vas  also  placed  on  this 
board  as  well  as  some  circuitry  to  supply  the  correct  docks  to  the  system. 

The  interface  between  the  modified  Volume  Surveillance  and  the  HP  hardware  consistad 
of  an  SBus  to  iEEE488  (QPIB  or  HPIB)  bus  converter  instafied  in  the  Sun.  This  converter 
allowed  for  direct  connection  of  the  Sun  and  the  built-in  HPIB  bus  of  the  HP  through  a  cable. 

The  HP  was  set  up  to  be  the  Controiler-ln-Charge  (CIC)  ai  the  bus  and.  therefore,  controlled 
the  data  transfer  from  the  Sun.  The  QPIB  bus  transfers  data  at  a  maximum  rate  of  1  Mbytefeec. 
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The  HP  still  received  digitized  and  other  data  from  the  NOVA.  Since  the  C30  only 
received  the  digitized  data,  it  could  only  do  a  portion  of  the  processing.  Therefore,  the  HP 
completed  the  remaining  processing.  After  the  C30  completed  its  processing,  it  passed  the  data 
to  be  displayed  and  other  data  necessary  to  complete  the  processing  to  tiie  HP. 

Software 

Although  the  original  Basic  code  was  translated  to  C,  some  modifications  were  made  for 
readability,  efficiency  and  to  correct  for  differences  between  the  two  languages.  Where 
possible,  variables,  algorithms  arfo  subroutine  names  were  kept  the  same  as  in  the  Basic  code. 
There  were  also  several  versions  of  the  software.  The  only  versions  covered  here  will  be  the 
original  translation  and  the  final  working  version. 

The  ori^nai  bandwidth  extrapolation  processing  translation  consisted  of  two 
subroutines.  The  first  subroutine  is  CaIcCoeff  vi^ich  calculates  the  complex  filter  coefficients. 
The  second  subroutine  is  Lnpred  which  linearly  extends  the  data  in  the  forward  and  backward 
direction  to  the  number  of  data  points  needed  in  the  extended  wavefonn.  The  FFT  code  used  was 
not  a  part  of  the  code  translate  from  Basic  to  C.  Two  subroutines  written  by  Jim  Evins  and 
Russell  Scott  in  C  performed  the  FFT  and  were  used  in  this  processing  as  well  as  the  stretch 
processing. 

The  processing  times  of  the  C  code  routines  were  found  using  the  C30  emulator.  The 
times  for  the  original  bandwidth  extrapolation  C  code  translation,  which  do  not  include  time  for 
data  input  and  output,  were  as  follows: 

CaIcCoeff  233  ms 

Lnpred  164  ms 

FFT  4.497  ms 


Total  Tune  401.5  ms 

The  processing  times  for  ttie  original  bandwidth  extrapolation  code  written  in  Basic  were 
approximately  5  seconds  and  3  seconds  for  the  subroutines  CaIcCoeff  and  LnPred  respectively. 
Therefore,  a  substantial  improvement  was  made  by  the  change  in  hardware  and  the  C  language 
conversion. 

The  original  strefoh  data  processing  translation  consisted  of  eight  subroutines.  The  list 
of  subroutines  along  with  a  description  of  what  the  subroutine  does  is  as  follows: 

1 .  Initialize  >  This  subroutine  initializes  some  constants,  sets  default  values  and 
calculates  the  phase  compensation  values. 

2.  Get_File  >  This  subroutine  gets  the  digitized  pulse  data,  performs  amplitude 
weighting  and  phase  compensation  on  the  data  and  calculates  the  values  of  other 
variables. 

3.  SignaLMag  -  This  subroutine  calculates  the  magnitude  of  two  pulses  of  digitized  data 
and  sends  this  data  out  to  be  displayed. 

4.  Rmotion.Calc  -  This  subroutine  performs  a  complex  multiply  between  the  first  and 
last  pulse  and  determines  the  range  offset 

5.  Stretch  >  This  subroutine  uses  the  compensated  data  from  GetLFile.  An  FFT  is 
performed  on  the  data  and  the  resulting  range  profiles  are  time  aligned  using  the 
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previously  determined  range  of^L  An  average  range  profile  is  found  by  magnitude 
averaging  the  range  profiles  of  each  pulse. 

6.  Range.Profile  -  This  subroutine  sends  the  logarithm  of  the  range  profiles  and  average 
profile  out  to  be  displayed. 

7.  Locate.Pmods  -  This  subroutine  calculates  a  threshold  used  to  determine  where 

propeller  modulation  occurs. 

8.  Decision  -  This  subroutine  is  called  within  Locate.Pmods  and  produces  a  target  class 
decision. 

The  processing  times  for  the  original  stretch  processing  C  code  translation,  which  do  not 
include  time  for  data  input  and  output,  were  as  foliows: 

Initialize  26.88  ms 

Get_File  45.56  ms 

Signal_Mag  4.852  ms 

Rmotion_Calc  119.9  ms 

Stretch  946.5  ms 

Range_Prof  33.85  ms 

Locate_Pmods  4.728  ms 


Total  Time  1.190  s 

The  processing  time  for  the  original  stretch  processing  code  written  in  Basic,  which  included 
the  eight  subroutines  listed  earlier  and  other  processing,  was  approximately  5  seconds. 

The  final  version  of  this  software  included  performing  the  bandwidth  extrapolation  on 
each  pulse  of  data  and  then  the  stretch  processing  on  the  group  of  extrapolated  puisM.  A  copy  of 
this  softwrare  appears  in  Appencfix  1.  To  meet  the  time  specifications,  the  processing  on  each 
S^up  of  pulses  must  be  completed  within  the  time  taken  for  one  scan  of  the  radar  which  is 
approximately  five  seconds.  For  the  case  of  seven  pulses,  the  time  taken  for  processing  is  equal 
to  (7  *  0.4015  s)  +  1.190  s  =  4.0005  s. 

As  mentioned  in  the  hardware  section  of  this  report,  the  C30  only  dki  a  portion  of  the 
processing.  Therefore,  some  subroutines  and  code  from  the  ori^nal  translation  are  omitted  in 
the  final  version  of  the  software.  Since  the  C30  does  not  receive  data  from  the  NOVA,  some 
needed  data  was  assumed  constant  by  the  C30.  The  data  that  was  assumed  constant  was  stored  in 
two  files,  one  for  all  the  variables  and  one  for  the  phase  compensation  data.  These  files  can  be 
erfited  on  the  Sun  computer  to  change  the  constants.  Any  changes  made  to  items  1  through  6  or 
the  number  of  phase  compensation  values  will  also  necessitate  changes  within  the  code.  The  data 
and  its  constant  value  was  as  follows: 


1.  Number  of  points  selected  from  full  waveform  333 

2.  Number  of  filter  coeffictente  150 

3.  Number  of  points  in  extended  waveform  1 000 

4.  Number  of  pulses  7 

5.  Number  of  samples  per  pulse  1 000 

6.  Number  of  points  in  FFT  4096 

7.  Sampling  rate  (MHz)  10 

8.  Bandwidth  200 

9.  Pulse  width  100 
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10.  Phase  compensation  data  (1000  floating  point  numbers) 

The  final  version  of  the  bandwidth  extrapolation  code  did  not  include  taking  the  FFT.  The 
final  version  of  the  stretch  processing  code  did  rK)t  include  any  of  the  code  relating  to  variables 
that  were  not  available  to  the  C30.  The  subroutines  Locate.Pmods  and  Decision  were  removed 
also.  Data  was  converted  to  16  bit  two’s  complement  integers  before  being  sent  to  the  displaying 
HP  as  a  stream  of  integers.  Floating  point  numbers  were  converted  to  integers  by  multiplying 
by  a  scaling  factor  and  truncating  the  fractional  portion  of  the  numbers.  The  data  sent  to  the  HP 
and  the  number  of  integers  it  was  composed  of,  in  parentheses,  appears  below  in  the  order  in 
which  it  was  sent: 

1.  Signal  magnitudes  of  two  pulses  (1000  per  pulse  =  2000)  -  scaling  factor  of  32. 

2.  Range  profile  of  each  pulse  (2000  per  pulse  s  14000)  -  scaling  factor  of  32768. 

(linear  data) 

3.  Average  range  profile  (2000)  -  scaling  factor  of  32768.  (linear  data) 

4.  Velocity  (1)  -  scaling  factor  of  100. 

5.  Magnitude  sum  of  one  pulse  (1)  -  scaling  toctor  of  100. 

6.  Range  oftoet  (variable  name  Pleat)  (1). 

7.  Peak  dwell  (1). 

Notes: 

-  Each  array  (1 ,  2  and  3)  are  sent  from  the  lowest  to  the  highest  index.  The  indices  for  the 
signal  magnitudes  are  1  to  1000.  The  indices  for  the  range  profiles  and  average  profile  are 
-999  to  1000. 

-  Arrays  marked  as  linear  data  are  not  ready  to  be  displayed  immediately;  the  logarithm  of 
the  array  is  displayed. 

Due  to  hardware  specifications  and  the  method  of  transmitting  data,  pairs  of  16  bit 
integers  were  packed  into  32  bit  words  in  the  C30.  A  byte  reversal  also  was  performed  on  the 
data  since  the  transputer  uses  ’litUe  endian”  byte  ordering  and  the  HP  uses  ”big  endian”  byte 
storage.  In  other  words,  the  transputer  stores  the  most  significant  byte  of  a  word  in  a  higher 
memory  address  than  the  least  significant  byte  whereas  the  HP  stores  the  most  significant  byte 
in  a  lower  memory  address.  During  transmission,  one  byte  at  a  time  was  sent  starting  vwth  the 
most  significant  byte.  Therefore,  the  receiving  system  waited  to  receive  two  bytes  for  each 
integer.  The  software  responsible  for  transmitting  the  data  across  the  IEEE  488  bus  was  run  on 
the  Sun.  When  all  of  the  cfota  was  received  by  the  Sun  and  ready  to  be  sent,  the  Service  Request 
(SRQ)  Pne  was  asserted.  The  HP  responded  by  allowing  the  Sun  to  send  the  data  over  the  bus. 

The  remaining  software  to  be  discussed  is  the  code  running  on  the  transputers  and  the 
main  program  to  control  the  processing  on  the  C30.  A  console  program  running  on  the 
transputer  resident  in  the  Sun  was  responsible  for  reading  the  data  from  tire  files  and  sending  it 
out  to  the  C30.  Several  transputers  along  the  network  were  responsible  for  passing  the  data  in 
the  correct  path  to  reach  the  transputer  on  the  SP  board  or  to  pass  the  data  to  the  GPIB  bus 
converter  for  transmission  to  tire  HP.  The  main  program  on  the  C30  was  responsible  for 
waiting  to  receive  the  data  read  from  the  data  files.  The  C30  polled  one  of  its  ports  for  an 
indication  ttiat  the  digitized  data  for  the  group  of  pulses  was  ready.  After  receiving  the  digitized 
data,  the  bandwidth  extrapolation  and  stretch  processing  occurred  and  the  necessary  data  was 
sent  out  The  process  started  again  with  the  C30  waiting  for  digitized  data  to  be  ready  and 
repeats  until  interrupted  externally. 
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Status 


Due  to  time  limitations,  this  real-time  processing  system  was  not  fully  implemented. 
The  software  was  fully  translated  from  the  Basic  language  to  the  C  language  and  was  tested  for 
correct  performance.  The  software  was  ported  to  the  final  system  located  at  the  Chesapeake  Bay 
Detachment  of  NRL  fitxn  the  testing  site  at  NRL  Washington.  The  modified  Volume  Surveillance 
hardware  successfully  received  digitized  data  from  the  Senrad  receiver  which  was  routed  to  the 
C30  processor.  The  C30  correctly  performed  the  processing  necessary  within  the  time 
specifications  and  sent  data  to  the  Sun  computer  to  be  relayed  to  the  HP  computer.  Some 
communication  between  the  Sun  and  HP  computer  across  the  IEEE  488  bus  was  implemented  but 
was  not  fully  functkmal.  The  only  remaining  task  is  to  modify  software  on  the  HP  computer 
which  would  correct  the  problems  with  the  interfsKse. 

Summary 

An  existing  system  using  NRL's  Senrad  radar  is  being  used  to  demonstrate  the  utility  of 
wideband  waveforms  in  radar  surveillance.  This  system  processes  data  using  a  stretch 
technique  to  yield  high  resolution  range  profiles  of  foe  target  To  provide  some  of  the  benefits  of 
greater  bandwidth  without  implementing  the  wideband  waveforms,  it  was  suggested  to  add 
bandwidth  extrapolation  processing.  The  addition  of  this  processing  necessitated  changes  to  the 
existing  system  for  real-time  processing.  The  changes  included  modifying  the  Volume 
Surveillance  system  hardware,  constructing  an  interface  between  the  existing  system  and  the 
Volume  Surveillance  processor  and  modifying  the  stretch  and  bandwidth  extrapolation  software 
for  a  new  processor.  The  changes  were  implemented  and  the  modified  hardware  operated 
correctly.  The  software  modifications  were  tested  successfully  with  recorded  arid  real-time 
data.  The  processing  times  met  the  time  specifications  necessary  for  real-time  operation. 
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^ '*«'***'***«*'*******  **************************•**************•****«**********'** 


Title: 

Stretch  Program 

File  Neune: 

str  .c 

* 

* 

Date: 

7/28/93 

it 

* 

* 

Purpose : 

To  do  stretch  processing  on  some  number 
1000  samples  each.  (See  block  diagram) 

of  pulses  of 

* 

* 

Target  HW: 

'mS320C30  Signal  Processor  Nodes 

* 

Author (s) : 

George  Linde 

* 

Athena  Caul  (rewrote  Basic  program  in  C) 

7/28/93 

it 

Revision: 

Test  code  reading  Chan  data  from  files. 

****************** ***************************************************«******y' 

/*  Note:  This  version  has  arrays  numbered  starting  from  0. 

To  be  used  for  creating  C30  code  */ 

flnclude  <stdlib.h> 
tinclude  <math.h> 
iinclude  •hw_sp.h* 
tinclude  ■dmsg.h* 
tinclude  ■msg.h* 
tinclude  “dp.h" 
tinclude  *dpprint.h* 
tinclude  'complex. mac* 

/*  includes  declaration:  struct  {  float  1,  q;  }  complex;  */ 
tinclude  'fft-h* 
tinclude  'xilinx.h* 
tinclude  ■ radar. h* 

/*****«***********  default  values  and  Constants  ************************/ 


tdef ine 

NPTS 

333 

/* 

t  of  points  selected  from  full  waveform*/ 

tdefine 

NCOEFF 

150 

/* 

t  of  filter  coefficients  */ 

tdef ine 

NEXT 

1000 

/* 

t  of  points  of  extended  array  */ 

tdefine 

NPTFFT 

4096 

/* 

t  of  points  in  FFT  */ 

tdefine 

NPULSES 

7 

/* 

Default  t  of  pulses  */ 

tdefine 

NSAMP 

1000 

/* 

t  of  samples  */ 

tdefine 

BW 

200 

/* 

signal  bandwidth  */ 

tdefine 

PW 

100 

r 

radar  pulse  width  */ 

tdefine 

SAMPRT 

1.E7 

/* 

sampling  rate  in  Hertz  */ 

tdefine 

FORWARD 

0 

/* 

Direct  Fourier  Transform  */ 

tdefine 

REVERSE 

1 

/* 

Inverse  Fourier  Transform  */ 

/************  global  variables  ***•********♦**•*****•**********************’ 

*  defined  global  for  use  of  external  functions  without  these  being  placed 

*  on  stack  (don't  need  to  pass,  initialized  to  zero) 


int  Chan_1181 (1000),  Chan_2l81 [1000] ; 
float  PhasetlOOO],  Spar[12]; 

float  SidOOO],  CoUOOO],  I_data[8]  [1000] ,  Q_data[8]  [1000]  ; 
int  M,  N,  Ns,  Np,  Pk_at,  Pk_dwell; 
float  P2,  Pv,  Mag  sum.  Vel; 
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y************************************************************************* 

*  CalcCoeff  calculates  the  conplex  filter  coefficients. 

*  Inputs  are:  NPTS,  NCOEPF,  Xit]. 
*************************************************************************/ 

static  void  CalcCoeff (int  Npts,  int  Ncoeff,  float  Pm[],  float  *Po, 
conqplex  X[],  complex  A[]) 

{ 

con«>lex  ’*Aa  =  {con5>lex  *)  malloc  (  550  *  sizeof  (complex)  );  /‘work  array*/ 

conplex  *B1  =  (complex  *)  malloc  (  550  *  sizeof (complex)  );  /*work  array*/ 

complex  *B2  =  (complex  *)  malloc  (  550  *  sizeof (complex)  );  /*work  array*/ 

complex  Xnom,  Temp,  Temp2; 
int  It,  M,  Mml,  Nmm,  Nml; 
float  Power,  Den; 


*Po  =  0.0; 

for  (It  =  0;  It  <  tpts;  It++) 

*Po  +*  (  CMAG2(  X[It]  )  ); 

*Po  /=  Npts; 

Itail  =  l^ts  -  1; 

B1[0]  =  XtO] ; 

B2[lbal  -  1]  =  X[Nmll; 

for  (It  =  1;  It  <  Nml;  It++)  { 

Bl(It]  =  X[It] ; 

B2[It  -  1]  =  Xlltl ; 

} 

for  (M  =  1;  M  <=  Ncoeff;  M++)  { 

Mml  =  M  -  1; 

Ntaim  s  l^ts  -  M; 
if  (  M  !=  1  )  { 

for  (It  =  1;  It  <=  Mml;  It++) 

Aa(It-l]  =  A(It-l] ; 
for  (It  =  1;  It  <=  Nnm;  It++)  { 

CONJ(  Aa(Hml-l],  Temp2  ); 

CMULT(  Temp2,  B2[It-l],  Temp  ); 

CSUB(  Bl[It-l],  Temp,  Bl[It  -  IJ  ); 

CMULT(  Aa[Mml-l],  Bltlt],  Temp  ); 

CSUB(  B2[It],  Temp,  B2[It-l]  ); 

) 

} 

Xnom.l  =  Xnom.q  =  Den  s  0.0; 
for  (It  =  1;  It  <=  Nora;  It++)  { 

CONJ(  Bl[It-l],  Temp2  ); 

CMDLT(  B2[It-l],  Temp2,  Temp  ); 

CADD(  Xnom,  Temp,  Xnom  ) ; 

Den  +=  (  CMAG2(BllIt-ll)  +  CMAG2 (B2llt-13 )  );  / *CADD (Den, Temp, Den) */ 

) 

if  (  Den  ==  0.0  ) 

A[lfinl}.i  =  A[Mml].q  »  O.O; 
else  ( 

CDIVS(  Xnom,  Den,  Temp  ); 

A[Mml].i  =  2.0  *  Temp.l;  /*  CMULT(  Two,  Temp,  A[Mml]  );  */ 

A[lftal]  .q  =  2.0  *  Temp.q; 

} 

Power  =  *Po; 
if  (  M  >  1  ) 

Power  =  Pm[M  -  2]; 

Pm[Mml]  =  Potrar  *  (  1.0  -  CHAS2(  A[Mml]  )  ); 
if  (  M  I*  1  ) 

for  (  It  =  1;  It  <=  Btal;  It++  )  { 

COHJ(  Aa[Mnl  -  It],  Temp2  ); 

CM0LT(  A[lfinl],  Temp2,  Temp  )  ; 

CSOB(  AaIIt-1],  Temp,  A[lt-1]  ) ; 
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} 

free 

free 

free 


Aa  ) 
(  B1  ) 
(  B2  ) 
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Lnpred  linearly  extends  the  complex  data  X  from  N1  points  to  N2 

*  points.  The  original  data  positioned  in  the  first  N1  elements  of 

*  array  X,  are  shifted  to  the  middle  of  X.  Forward  and  backward 

*  predictions  are  done  until  the  total  #  of  data  points  is  N2. 

*  Note:  the  N2  points  contain  the  original  N1  points. 

A***********)^******************************** **************************** ^ 

static  void  Lnpreddnt  Nl,  int  N2,  int  Ncoeff,  complex  X[],  coiif>lex  A[]) 

{ 

int  I,  J,  K,  LI,  L2.  N3; 
con^lex  Temp,  Ten%>2; 

LI  =  N2/2  -  Nl/2; 

L2  =  N2/2  +  Nl/2; 
if  (  (N1  %  2)  ==  1  )  { 

LI  +=  1; 

L2  +=  1; 

} 

for  (I  =  0;  I  <  Nl;  I++  )  { 

J  =  N1  -  I  -  1; 

K  =  L2  -  I  -  1; 

XtK]  =  X[J] ; 


N3  =  N2  -  L2;  /*  Do  for%irard  predictir*,  ’/ 

for  (I  =  0;  I  <  N3;  I++  )  { 

J  =  L2  +  I; 

X[J] .i  =  XtJ] .q  =  0.0; 

for  (K  =  0;  K  <  Ncoeff;  K++  )  { 

CMULTi  A[K],  X[J-K-1],  Ten?)  ); 

CADD(  XtJ],  Ten?),  XtJ]  ); 

} 

} 

for  (I  =  0;  I  <  LI;  I++  )  {  /*  Do  backward  prediction  */ 

J  =  LI  -  I  -  1; 

XtJ] .i  =  XtJ] .q  =  0.0; 

for  (K  =  0;  K  <  Ncoeff;  K++  )  { 

CONJt  AtK],  Ten?)2  ); 

CMULTt  Ten?)2,  XtJ+K-fl],  Temp  ); 

CADD(  XtJ],  Ten?),  X[J]  ); 

} 

} 

} 


/*  set  up  limits  for  loops  */ 


/*  shift  original  data  to  middle  of  Xt]  */ 
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^  *******************  *******«******'*'*«««««*'*««**************************•«* 

*  Initialize 

*************************************************************************/ 

void  Initialize (  ) 

{ 

int  I,  J; 

float  K,  SI,  Ys,  Yl,  Phr,  G(4]; 
double  loglOO,  sin(),  cob(); 


Np  =  NPULSES;  /* 
Ns  =  NSAMP;  /* 
SI  =  BW  /  PW;  /* 
P2  =  (SAMPRT  /  SI)  /2033.4;  /* 
M  =  (int)  (  loglO( (float)  NPTPFT) 
N  =  1  «  M;  /* 


Nunber  of  pulses  */ 

Nutaber  of  sasoples  */ 

Chirp  slope  */ 

Scale  in  feet  */ 

/  logl0(2.0)  +  0.5  ); 

N  =  points  in  FPT  =  2''H  */ 


G[l]  s  0.29265;  /*  phase  conpensation  */ 

G[2]  =  -0.0157838; 

G[3]  s  0.00218104; 

for  (J  =  0;  J  <  Ns;  J++  )  (  /*  read  in  baclcwards  */ 

K  =  ( (float) J  -  ((float)  N8)/2.)  /  ((float)  Ns) ; 

Ys  =  0.0; 

Ys  =  Gtl]  *  cos(K*2*M_PI)  +  Gt2]  *  co8{K*4*M_PI)  +  G[3]  *  cos (K*6*H_PI) ; 
Yl  *  0.6  +  1.2  *  Ys; 

Phr  =  (-  Phase (J])  *  K_PI  /  720.0;  /*  {-tp)/4  in  radians  */ 

SitJ]  =  sin(Phr)  *  Yl; 

Co[J]  =  cos (Phr)  ♦  Yl; 

} 

} 


/************************************************************************* 

*  Get_File  gets  the  stretch  data  and  performs  lunplltude  weighting 

*  and  phase  conpensation  (using  values  calculated  in  Initialize) . 
*********«*************•********************************************«****/ 

void  Get_Pile(  ) 

{ 

int  I,  J,  tpl,  tp2; 
float  S,  SI,  Bw,  Pw; 
double  sqrt ( ) ; 


Np  =  Spar [3] ; 
Ns  =  Spar [4] ; 
N  s  Spar[5] ; 


S  s  Spar [6] ; 

/• 

sanpling  rate  in  MHz  */ 

Bw  =  Spar[7] ; 

Pw  s  Spar [8] ; 

P2  =  S  *  1000000; 

/* 

sanpllng  rate  in  Hz  */ 

SI  =  Bw  /  Pw; 

/* 

chirp  slope  */ 

P2  =  (P2  /  SI)  /  2033.4; 

/*  scale  in  feet  */ 

for  (  I  =  0;  I  <  Np;  I++  )  { 

/* 

read  all  pulses  in  dwell  */ 

for  (J  =  0;  J  <  Ns;  J++  )  { 

/* 

read  in  baclcwards  to  correct 

tpl  =  Chan_l[I]lJ]  /  28.; 
tp2  =  Chan_2[I](J]  /  28.; 

Chan^llUlJ]  =  tpl; 

Chan^2II][J]  =  tp2; 

I_data[I][J]  =  (  ((float)tpl)  *  S1[J1  +  Co[J]  *  ((float)tp2)  )/1000.0; 
Q_datalIl[J]  =  {  {(float)tp2)  *  SiU]  -  CoUl  *  ((£loat)tpl)  )/1000.0; 

) 

} 

} 
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*  Signal_Hag  calculatss  tha  signal  magnitude  of  the  I  Q 

*  and  plots  it  to  a  file. 


void  Signal_Mag(  dmsg  ) 

CMsg  *diiisg; 

{ 

int  I,  Yt,  Yp,  Y,  cnt; 

float  *Magl  =  (float  *)  malloc  (  1000  *  sizeof (float)  )j 
float  *Mag2  =  (float  *)  malloc  (  1000  *  sizeof (float)  ) ; 
double  sqrt ( ) ; 

for  (I  =  0;  I  <  Ns;  I++  )  { 

Yt  =  Chan_2[3] (11; 

Yp  =  Chan_l(3] [I]; 

Magi [I]  =  (float)  (sqrt ( (float)  (Yt*Yt  +  Yp*Yp) )  *  1.41); 
Yt  *  Chan_2(Np-l] [I] ; 

Yp  =  Chan_l[Np-l] [I]; 

Mag2[I]  =  (float)  (sqrt ( (float)  (Yt*Yt  +  Yp*Yp) )  *  1.41); 
if  (  I  <  20  ) 

DP_PrintfC  Magl[%d]  =  %d\n*,  I,  (int)  (Magl[I]  *  32); 


cnt  =  0; 

for  (  I  =  0;  I  <  Ns;  I  =  I  +  2  )  { 

Y  *  (int)  (Magi [I]  *  32); 

Yt  =  (  (Y  «  8)  &  OxFFOO  )  t  (  (Y  »  8)  &  OxOOFF  ) ; 

Y  *  (int)  (Magl[i+l]  •  32); 

Yp  =  (  (Y  «  8)  &  OxFFOO  )  |  (  (Y  »  8)  fc  OxOOFF  ) ; 

^  dinBg->8igmag.inaglcnt++l  =  (  (Yp  «  16)  |  (Yt  6  OxFFFF)  ); 

for  (  I  =  0;  I  <  Ns;  I  =  I  +  2  )  { 

Y  =  (int)  (Mag2tll  *  32); 

Yt  =  (  (Y  «  8)  Sc  OxFFOO  )  |  (  (Y  »  8)  A  OxOOFF  )  ; 

Y  =  (int)  (Mng2[I+l]  *  32); 

Yp  =  (  (Y  «  8)  A  OxFFOO  )  |  (  (Y  »  8)  A  OxOOFF  ); 
dinsg->sigmag.mag[cnt+'t-]  =  (  (Yp  «  16)  I  (Yt  A  OxFFFF)  ); 
if  (  I  >  996  ) 

DP_Printf ("sigmagtAdl^Ad  sigmag[%d]s%d\n*, I, (int) (Mag2[Z]*32) , 
^  I+l, (int) (Mag2[I+ll  *  32)); 

DP_WriteOut(  DMSG_TYPEL.SIGMAG_PATA, 

sizeof (EMsg_SigmagData) ,  (Msg  *)dmsg  ); 

Mag_8Uin  s  0.0; 

for  (I  =  0;  I  <  Ns;  I++  ) 

Mag_sum  Magl[I]; 

Mag_8um  /=  ((float)  Ns); 
free (  Magi  ) ; 
free (  Mag2  ) ; 
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/**»***************************************************************** 

*  Riiiotion_Calc  finds  the  product  of  t%io  pulses  and  determines  the 

*  range  offset. 


void  Bin0tion_Calc  (  ) 

{ 

int  I,  J,  M,  Lmi; 

/*  Subscript  of  Yr  is  actually  -128  to  127  so  subtract  128  */ 
float  Yoffset,  Pvf,  XI,  X2,  Yl,  Y2,  Tim; 
float  *Yr  =  (float  *)  malloc  (  257  *  sizeof (float)  ); 
complex  *tbl  =  (complex  *)  malloc  (  2050  *  sizeof (complex)  ); 

/*  Chan  is  R_float  and  I_float  */ 
complex  *Chem  =  (complex  *)  malloc  (  4096  *  sizeof (complex)  ); 
double  loglO ( ) ,  sqrt ( ) ; 

for  (  I  =  0;  I  <  4096;  I++  ) 

Cheui  [  I  ]  .  i  =  Chan  [  I  ]  .  q  =  0.0; 
for  (  I  *  3;  I  <  997;  I++  )  { 

XI  »  ((float)Chan_l(0] (ID  /  1000.; 

X2  s  ((float)Chan_l(Np-l] (ID  /  1000.; 

Yl  s  ((float)Chan_2(0] (ID  /  1000.; 

Y2  *  ((float)Chan_2(Nt>-l]  (ID  /  (-1000.); 

Chandl.i  =  X1*X2  -  Y1*Y2; 

ChandJ.q  =  X2*Y1  +  X1*Y2; 

) 

M  =  12;  /*  N  =  2''M  points  in  the  FFT  =  4096  */ 

I  =  1  «  M; 

fft_init(  I,  FORWARD,  tbl  );  /*  FFT_JPlonting ( ) ;  */ 

/*  complex  array  Chan  is  used  for  both  input  and  output  */ 
fft2(  Chan,  I,  FORWARD,  tbl  ); 


/ 


for  (  I  =  128;  I  <=  255;  I++  )  /•  from  0  to  127  */ 

Yrd]  -  20.0  *  ((float)  logl0(  CMAG (Chan (1-128])  l.OE-15  )); 

for  (1=0;  I  <=  127;  I++  )  /*  from  -128  to  -1  */ 

Yrd]  =  20.0  *  ((float)  logl0(  CHAO  (Chan  (3968+1  ] )  +  l.OE-15  )  ); 

ijni  =  0;  /*  find  maximum  in  Yr(]  */ 

Yoffset  =  Yr[0]; 

for  (  I  =  1;  I  <=  255;  I++  )  /*  from  -128  to  127  */ 

if  (  Yrd]  >  Yoffset  )  { 

Yoffset  =  Yrd]; 

Lmi  =  I; 

) 

Tim  =  (Np  -  1)  *  0.0031284; 
if  (  (Lmi  >  255)  |  |  (Imii  <  2)  ) 

Lmi  =  128; 

XI  =  Yr(Imi-l]  -  Yoffset; 

X2  =  Yr(lmi+1]  -  Yoffset; 

Yl  =  Yrdrni]  -  Yoffset; 

Pv  =  ((float) Lmi) -128.  +  0.5  *  (XI  -  X2)  /  (XI  +  X2  -  2*Y1); 

Pvf  =  Pv  *  P2  /  N; 

Vel  =  Pvf  /  Tim  *  3600.  /  6076.1155;  /♦  round  this  to  two  digits  */ 

free (  Yr  ) ; 
free (  tbl  ) ; 
free (  Chan  ) ; 
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^*******'*'***************************************************************** 

*  Stretch  uses  the  compensated  data  from  Get_file.  Stretch  performs  an 

*  FFT,  time  alignment  using  values  from  Rmotloi}_Calc  and  does  a 

*  magnitude  average  on  the  range  profile. 
******«********************************«*«***••********«**********'*******/ 

void  Stretch(  R_profile,  Profile_avg  ) 
float  R^rofile[]  [2401] ,  Proflle_avg[ ] ; 

{ 

int  1,  J,  L,  M,  Last.dwell,  Vel_offset,  T^offset; 
float  max,  Dwell_mag,  sum,  X; 

float  *p_tmp  s  (float  *)  malloc  (  4096  *  slzeof (float)  ) ; 
complex  *tbl  =  (complex  *)  malloc  (  2050  *  sizeof (complex)  ); 

/*  data  is  Refloat  and  I_float  */ 

complex  *data  =  (coiiv>lex  *)  malloc  (  4096  *  sizeof  (complex)  ); 
double  sqrt ( ) ,  loglO ( ) ; 

Last_dwell  =  0; 
for  (  I  =  0;  I  <  Np;  I++  )  { 
for  (  J  =  0;  J  <  Ns;  J++  )  { 
data[J].i  =  Q_data[I] [J] ; 
data[J].q  =  I_data[I] [J] ; 

} 

for  (  J  =  Ns;  J  <  4096;  J++  ) 
data[J].i  =  data[J].q  =  0.0; 

M  =  12;  /*  N  *  2'‘M  points  in  the  FFT  =  4096  */ 

L  =  1  «  M; 

fft_init(  L,  FORWARD,  tbl  ) ;  /*  FFT^loating  ( ) ;  */ 

/*  complex  array  data  is  used  for  both  input  and  output  */ 
fft2(  data,  L,  FORWARD,  tbl  ); 

/*  data[0-4095]  =  ICtmp(2048-4095, 0-2047]  originally  [1-2048,-2047-0]  */ 
for  (  J  =  0;  J  <  2048;  J++  ) 

R_tmp[J+2048]  -  CaiAG(  data[J]  );  /*  from  2048  to  4095  */ 

for  (  J  =  2048;  J  <  4096;  J+-^  ) 

R_tmp[J-2048]  =  CMAG(  data[J]  );  /*  from  0  to  2047  */ 

if  (  I  ==  0  )  {  /*  locate  pealc  6  velocity  offset  */ 

X  =  0; 

max  =  R_tn«>[0];  /*  locate  max  R^tmp  value  */ 

for  (  J  =  1;  J  <=  4095;  J++  )  /*  from  -2047  to  2048  */ 
if  (  R^tmp[J]  >  max  )  { 
nmx  =  R^tmp[J] ; 

X  =  J; 

} 

P)C_at  =  X  -  2047; 

} 

Vel_offset  =  (int)  (0.5  +  Pv  /  6.  *  (I+D); 

T_off8et  =  0; 

T_offset  =  Vel_offset  +  Plt_at; 
if  (  I  ==  6  ) 

DP_Printf('-  Str  -  Voff  %d  P)«L.at  %d\n*,Vel_off8et,  PK_at)  ; 
for  (  J  =  0;  J  <=  2400;  J++  ) 

Rj>rofile[Il  [J]  =  IL.tmp[847+T_of fset+J] ;  /*  offset  data  curray  */ 

Dwelljmag  =  0.0; 
for  (  J  =  0;  J  <=  4095;  J++  ) 

Dwell_mag  +-  R_tmp[J]; 

if  (  Dwelljnag  >  Last.dwell  )  {  /♦  save  dwell  *  of  peak  dwell  to  */ 
Last.dwell  =  Dwell_mag;  /*  use  in  res  measurement  */ 

Pk^dwell  =  I; 

1 


) 
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free  (  I<_tinp  )  ; 
free(  tbl  ) ; 
free (  data  ) ; 

/*  averaged  -  euni  total  of  profiles  divided  number  of  pulses  */ 
for  (  J  =  0;  J  <=  2400;  J++  )  { 
sum  =  0.0; 

for  (  1=0;  I  <  Np;  I++  ) 
sum  +=  Rj)rofile[I)  [J]  ; 

Profile_avg[J]  =  sum  /  ((float)  Np) ; 
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/* 

* 

* 


Range_Proflle  plots  the  range  profile  for  each  pulse  separately 
emd  the  profile  average  of  all  the  pulses. 


void  Range_Profile(  R_profile,  Profile_avg,  dmsg  ) 
float  R_profile[] [2401],  Profile_avg( ] ; 

DMsg  *dinsg; 

{ 


/ 


int  I,  J,  a,b,c,d; 

float  X,  Y; 

double  loglOO; 

DP_Printf( 'PEAK  LEVEL  :  %6. If  dB\n', (float)  (20.*logl0(Profilo_avg[1200] ) ) ) ; 

DP_Printf  ( -RANCSE  OFFSET  :  %5d  FT\n-.  (int)  (( (float) PK_at)  *  0.6)  ); 

d  s  0; 

for  (  1  =  0;  I  <  I++  ) 

for  (  J  =  -999;  J  <=  1000;  J  =  J  +  2  )  { 
c  =  (int)  (  R_profile[I] (J+12001  *  32768); 
a  =  (  (c  «  8)  &  OxFFOO  )  |  (  (c  »  8)  6  OxOOFF  ); 
c  =  (int)  (  R_profileII] (J+1201]  *  32768); 
b  =  (  (c  «  8)  &  OxFFOO  )  I  (  (c  »  8)  fc  OxOOFF  ); 
dinsg->rgprof  .rgld++]  =  (  (b  «  16)  |  (a  fc  OxFFFF)  ); 

} 

DP_WriteOut(  DMSG_TYPE_RGPROF_EATA, 

sizeof (DM8g_RgprofData) ,  (Msg  *)dm8g  ); 


d  =  0; 

for  (  J  =  -999;  J  <=  1000;  J  =  J  +  2  )  { 
c  =  (int)  (  Profile_avg[j4-1200]  *  32768); 
a  =  (  (c  «  8)  &  OxFFOO  )  |  (  (c  »  8)  fc  OxOOFF  ); 
c  =  (int)  (  Profile_avglJ+12011  ♦  32768); 
b  =  (  (c  «  8)  &  OxFFOO  )  I  (  (c  »  8)  &  OxOOFF  ); 
dinsg->avgprofpl.avg[d++]  =  (  (b  «  16)  |  (a  &  OxFFFF)  ); 

} 

} 


/**************«****«**««*•*«*«*«•«*««*««««**««*•**•*«***«• 

*  c30_init  is  copied  from  /home/volume/src/c30/main.c 

*  This  function  removes  wait  states  and  enables  the  cache 


static  void  c30_init() 

{ 

*PRIMARY_CTL  t=  OxOOOOlflf 
*PRIMARY_CTL  &=  OxOOOOOOff 
*EXPANSION_EUS  6=  OXOOOOOOlf 
*ADDIL.S1_XMIT_CNTL_RBG  »  0x2 
asmC  or  800h,  st*  ); 

) 


/*  force  0  wait  states  P.Bus 
/*  No  ban]c  switching  */ 

/*  force  0  wait  states  E.Bus 
/*  Disable  Ext  Rom,  enable  Ram  */ 

/*  enable  cache  */ 


/ 


*/ 

*/ 
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malnO 

{ 

static  con^lex  Xi[1250]; 

static  float  Pm[550];  /‘array  of  updated  error  power*/ 

static  float  Po;  /‘real  variemce  of  the  data*/ 

static  conplex  A[550];  /‘array  containing  filter  coefficients*/ 

static  int  J,  Ml,  Mcoeff,  M2,  Of,a,b,c; 

static  int  I,  ctl,  ct2,  ct3,  count,  cnt; 

static  int  k,  1,  i_sanple,  q_aainple,  ‘intptr; 

static  float  Prof ile_avg[2401] ,  Rjprof ile[9] [2401] ; 

static  Radar_CaptureBuf fer  ‘rcb; 

IXlsg  ‘dmsg  =  (DMsg  *)  inalloc(  8lzeof(DMsg)  ); 

Msg_Length  length; 

M8g_Type  type; 

c30_init() ; 

Xilinx_Inlt() ;  /*  Initialize  Xlllnx  interface  */ 

I  =  Ct2  =  count  =  0; 
ctl  =  ct3  =  999; 

*ADDR_LED_PORT  =  -0x03; 

for(;;)  { 

switch  (  type  =  DP_ReadIn(  &length,  (Msg  ‘)dinsg  )  )  { 

case  DMSG_TYPE_XILIMX: 

handle_xilinx  (  (aisg_Xilinx  *)  dmsg  ); 

rcb  =  Radar_Init() ;  /*  Initiliaze  DMA  interface  */ 

DP_Prlntf(  'SP:  INFO:  Xilinx  program  received  emd  handled. \n*  ); 
break; 

case  E»!SG_'m»E_PHASEJ)ATA; 

Phase[ctl — ]  =  dm8g->phdata.p;  /*  read  in  backwards  temp  */ 
*ADDR^LED_PORT  *  -(++count); 
if  (  ctl  <  0  ) 
break; 

case  EMSG_TYPE_SPAR_11ATA: 

Spar  [ct2-:-+]  =  dmsg->8pdata.8p; 

*ADDR_I.ED_PORT  *  -(count++); 
if  (ct2  ==  9)  { 

Np  =  Spar [3] ; 

Ms  =  Spar [4] ; 

DP_Printf (‘All  of  the  Spar  data  receivedXn*) ; 

} 

break; 

case  0MS6_TXPE_CHAM_DATA: 

Chan_l[I] [ct3]  -  dmsg->chdata.one;  /*  read  in  backwards  temp  */ 
Cheuk_2  [  I  ]  [  ct3 — ]  =  dmsg->chdata .  two ; 

*ADDR_LED_PORT  =  -(count++); 
if  (  (ct3  <  0)  &&  (I  <  I^)  )  ( 

I++; 

ct3  =  999; 

} 

if  (  (ctl  <  0)  &&  (ct2  »  9)  &&  (I  =»  Mp)  )  { 

Ml  =  Spar[0];  /*  Pick  middle  Ml  points  */ 

Mcoeff  s  Spar[l]; 

M2  -  Spar [2];  /*  f  of  points  of  extended  array  */ 

DP_Printf (‘The  nuniber  of  data  points  is  %d\n‘,MPTS); 

DP_Prlntf ( ‘The  nuniber  of  coefficients  is  %d\n\n*, MCOEFF) ; 
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} 


for  (  cnt  =  0;  cnt  <  1;  cnt++  )  {  /*  infinite  loop  eventually  */ 

for  (  I  =  0;  I  <  Nt>;  I++  )  { 

Po  =  0.0; 

Of  =  N1  /  2; 

if  (  (N1  %  2)  *=  1  ) 

Of  -*  2;  /*  To  match  original  mveform  location  */ 

for  (  J  *  0;  J  <  Nl;  J++  )  { 

Xi(J].i  =  (float)  (Chanel [I] I  N2/2  -  Of  +  J  ]  ); 

Xi(J].q  =  (float)  (Chan_2tllt  N2/2  -  Of  +  J  )  ); 

> 

CalcCoeff (NPTS,  NCOEFF,  Pm,  APo.  Xi,  A) ; 

Lnpred(  Nl,  N2,  NCOEFF,  Xi,  A  ) ; 
for  (  J  =  0;  J  <  Ns;  J++  )  { 

Chan_l[Il[J]  =  (int)  (XiUl.l); 

Chan_2[I](J]  =  (int)  (Xi(Jl.a); 

} 

} 

Initialize (  ): 

OetJFile(  ); 

Signal_^g(  dmag  ) ; 

Rmotion_Calc (  ) ; 

Stretch(  R_profile,  Profile_avg  ); 

Range_Profile(  ILprofile,  Profila_avg,  dmsg  ); 

c  =  (int)  (Vel  *  100); 

a  =  (  (c  «  8)  &  OxFFOO  )  |  (  (c  »  8)  A  OxOOFF  ); 

c  s  (int)  (  Mag  sum  *  100  ); 

b  =  (  (c  «  8)  6  OxPFOO  )  I  (  (c  »  8)  6  OxOOFF  ); 

clmsg->avgprofpl . velmag  =  (  (b  «  16)  |  (a  &  OxFFFF)  ) ; 

a  =  (  (Pk^at  «  8)  &  OxFFOO  )  |  (  (P)t_at  »  8)  A  OxOOFF  )  ; 
b  =  (  (Pl«L.dwell  «  8)  A  OxFFOO  )  1  (  (P)t^dwell  »  8)  A  OxOOFF  ); 
dm8g->avgprof pi. peaks  =  (  (b  «  16)  )  (a  A  OxFFFF)  ); 
*ADDILJiED_PORT  =  -0x55; 

DP_WriteOut  (  DMSa_TYPELjAVGPROF_PItJaATA, 

sizeof  (OHsgL_AvgprofFlusData) ,  (Msg  *)dmsg  ); 

DP_Printf CXnVel  =  %8.2f  as  int  %d\n',  Vel,  (int) (Vel*100)  ); 
DP_Printf ( ■Mag_sum  =  %8.2f  as  int  %d\n*,Mag_sum, 

(int) (Mag_s\m  *  100)  ); 

DP_Printf  ('PjcL-at  =  %d\n",  Pk_at  ); 

DP_Printf CPk^dwell  *  %d\n-,  Pk_dwell  ) ; 

} 

} 

break; 

default: 

DP_Printf(*  ERROR:  Bad  message  type  received  (t*%02x  l=%d).\n*, 
type,  length  ) ; 

break; 

} 


static  void  handle_xilinx(  EMsg  Xllinx  *xilink_ptr  ) 

{ 

int  i; 

/*  Start  configuration  xfr  mode  *! 
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DP_Prlntf CSP;  ERROR:  Failed  xilinx  start. \n- ); 

Led_Print (  OxFl  ) ; 

/*  Send  configuration  data  */ 
for  (  i  =  0;  i  <  xilinx_ptr->n;  i++  ) 
Xilinx^Write(  xilinxj>tr->wordti]  ); 

Led_Print {  0xF2  ) ; 

/*  Make  sure  xilinx  chip  is  fully  configured  */ 
if  (  !Xllinx_Done()  ) 

DP_Printf CSP:  ERROR:  Failed  xilinx  done.\n-); 
t*ed_Print  (  OxFO  )  ; 
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